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Abstract 



We study six-point gluon scattering amplitudes in J\f = 4 super Yang-Mills the- 



ory at strong coupling by investigating the thermodynamic Bethe ansatz equations 
of the underlying Z^-symmetric integrable model both analytically and numerically. 
By the conformal field theory (CFT) perturbation, we compute the free energy part 
of the remainder function with generic chemical potential near the CFT/small mass 
limit. Combining this with the expansion of the Y- functions, we obtain the remain- 
der function near the small mass limit up to a function of the chemical potential, 
which can be evaluated numerically. We also find the leading corrections to the re- 
mainder function near the large mass limit. We confirm that these results are in 
good agreement with numerical computations. 
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1. Introduction 

Recently there has been much progress in computing the area of minimal surfaces 
in AdS space, whose boundary is made of light-like segments [IH7J. The surface 
area corresponds to the expectation value of the Wilson loop along the same contour 
and is dual to the gluon scattering amplitude in planar M = 4 super Yang-Mills 
theory at strong coupling. The minimal surface area has deviation from the Bern- 
Dixon-Smirnov (BDS) conjecture [8] by an amount of the remainder function, which 
is shown to exist for the n(> 6)-point amplitudes J9j[10]. Determination of the 
remainder function is a key step toward establishing the correct analytic formula of 
the gluon scattering amplitude. 

In [3J, the equations for determining the minimal surface in AdSr, space are shown 
to be the SU(4) Hitchin equations with certain constraints and boundary conditions. 
The area is characterized by the Stokes data of the asymptotic solutions of the 
associated linear system. The Stokes data obey certain functional equations and the 
minimal surface area is evaluated by solving the integral equations associated with 
them. For the 6-point gluon scattering amplitudes, the functional equations and the 
integral equations turn out to be the Y-system [11] and the thermodynamic Bethe 
ansatz (TBA) equations [12] of the Z4-symmetric (or A3) integrable model [13—15], 
respectively. The minimal surface area is evaluated by the free energy of the model. 

This solution has been generalized to the minimal surface with an n-sided light- 
like polygonal boundary in AdS5 [5] and in AdS3 [5j[6]. In particular, Alday, Malda- 
cena, Sever and Vieira proposed the Y-system and the TBA equations for the n-sided 
polygonal solution and expressed the area in terms of the TBA system. In [6], the 
present authors noticed that the TBA system is that of the homogeneous sine- Gordon 
model [To] . 

The free energy for a two-dimensional integrable model with purely elastic S- 
matrix is obtained by solving the TBA system. It is very difficult to solve the TBA 
system exactly but one can investigate its solutions in two limits. In the UV (or high 
temperature or small mass) limit, it is described by a certain conformal field theory 
(CFT). In the IR (or low temperature or large mass) limit, on the other hand, it 
becomes a system of a free massive theory. In the UV limit, the free energy turns 
out to be the central charge of the CFT. In the case of the homogeneous sine-Gordon 
model, the relevant CFTs are the generalized parafermion CFTs [17] . A particularly 
interesting feature of the TBA system is that it contains chemical potential which 
arises from the monodromy of the asymptotic solution of the linear system. Its 



solution with chemical potential near the CFT point has not been well studied. 

While analysis in the above described limits reveals several principal characters 
of the amplitude, evaluation of the remainder function apart from the limits is also 
of great significance in studying how the amplitude depends on the gluon momenta. 
The purpose of this paper is to analyze such momenta dependence of the remainder 
function. To do this, we investigate the TBA system perturbatively both near the 
UV and IR limits. In the present work we will focus on the 6-point amplitude 
for simplicity, where the relevant CFT is the Z^-parafermion CFT. The solution of 
the TBA system without chemical potential around the CFT point was studied by 
Klassen and Melzer [IB] . In this paper we will investigate the perturbative solution 
around the CFT point including chemical potential. We also study the TBA system 
near the IR limit and present the first correction to the remainder function. 

Our results provide an analytic form of the remainder function for the 6-point 
scattering amplitude at strong coupling away from the CFT and the IR point. 
This paper also demonstrates that the unexpected connection between the four- 
dimensional super Yang-Mills theory and the two-dimensional integrable models 
discovered in [21E] enables one to compute the amplitude at strong coupling in two- 
dimensional approaches. From a point of view of the study of integrable models, our 
discussion provides a concrete example of an analysis of TBA systems with chemical 
potential. 

This paper is organized as follows. In sect. 2, we review the construction of the 
Hitchin equations of the minimal surface with a 6-sided polygonal boundary in AdSs 
and the related Y-system and TBA equations. In sect. 3, we study the CFT limit of 
the TBA system and evaluate the free energy perturbatively. In sect. 4, we study the 
remainder function around the CFT limit. In sect. 5, we examine the large mass limit 
and obtain the correction to the remainder function. We conclude with a discussion 
in sect. 6. 

2. Review of TBA system for six-point amplitudes 

In this section we review the TBA system for the six-point gluon scattering ampli- 
tudes. We basically follow the references [HJIS]- 

2.1. Classical string solutions with a null polygonal boundary 

Alday and Maldacena proposed a method of computing gluon scattering amplitudes 
in M = 4 super Yang-Mills using AdS/CFT correspondence [lj. According to their 



proposal, scalar magnitude of n-gluon MHV scattering amplitudes can be evaluated 
in the strong coupling limit by computing the area of corresponding classical open 
string solutions. The string solutions are minimal surfaces whose boundary is a 
polygon located on the boundary of AdS$. The polygon consists of n null edges 
given by the n momenta of incoming gluons. 

To be concrete, classical string solutions under consideration are realized as a map 
from their world-sheet to AdS^. We consider Euclidean world-sheet parametrized by 
z, z. In terms of the global coordinates, AdS$ is expressed as a hypersurface 

X-X = -(X- 1 ) 2 - (X ) 2 + (X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 + (X 4 ) 2 = -1 (2.1) 

in IR 2 ' 4 . Classical string solutions X(z, z) satisfy the equations of motion 

ddX - (dX ■ dX)X = (2.2) 

and the Virasoro constraints 

OX ■ dX = dX ■ dX = 0. (2.3) 

The boundary condition is expressed in terms of the Poincare coordinates (r, x M ) 
given by 

X» = — , /i = 0,l,2,3, (2.4) 

r 

X- 4 + X 4 = - X' 1 -X 4 = +X % r (2.5) 

The solutions which correspond to n-gluon amplitudes end on an n-sided polygon at 
the AdS boundary r = 0. The vertices x%, . . . , x n of the polygon are separated by 
the gluon momenta ki,...,k n as 

X J ~ x i + i = k r (2-6) 

This type of boundary condition is neatly characterized in terms of the generalized 
sinh-Gordon potential [Tj5]. Let us introduce the following notation 



P{z) = d 2 X ■ d 2 X, P(z) = B 2 X ■ B 2 X, (2.8) 

a(z,z) = a(z,z) - i log P(z)P(z), (2.9) 

where a, a are real and P(z) is shown to be analytic in z. For the simplest four-cusp 
solution, P(z) = 1 and at = on the whole z-plane. For ra-cusp solutions, P(z) is 



a polynomial of degree n — 4 and a — > for \z\ — > oo. The latter condition reflects 
the fact that n-cusp solutions have the same asymptotics with the four-cusp solution 
around each cusp. 

In this paper we will concentrate on the case of n = 6, where P(z) is quadratic. 
One can choose a gauge 

P(z) = z 2 -U } U EC (2.10) 

by a suitable redefinition of the world-sheet coordinate. 

2.2. Hitchin system 

The equations of motion (12. 2 p and the Virasoro constraints (12.31) can be rephrased 
as SU(4) Hitchin equations [3J. To see this, let us consider a moving-frame basis 
spanned by go — X, g*4 = e~ a ^ 2 dX, g*5 = e~ a ^ 2 dX and the other three complementary 
orthonormal vectors gj, I = 1,2,3. The evolution of the basis q = (go, • • • , qV) is 
described by a set of linear differential equations 

dq = -A z q, dq = -A- Z q. (2.11) 

One can decompose the connection into two parts A = A + $, in such a way that A 
rotates (go, gi,g*2, gl?) and (g4,g*5) separately among themselves while $ mixes them. 
The flatness condition of (12.111) then takes the form of the Hitchin equations 

D z $- Z = 0, D- Z § z = 0, (2.12) 

[D z ,D s ] + [$ z ,®,} = 0, (2.13) 

where D z = d + [A z , ], D z = d + [A s , ]. These equations are equivalent to the 
equations of motion (12. 2p and the Virasoro constraints (12. 3p . One can write down 
the same equations in the spinor basis, where A Z ,A Z , $ 2 , <3> f now represent 4x4 
matrices. We do not need their explicit form [5] here, but an important fact is that 
these A and $ corresponding to string solutions satisfy additional constraints 

CA^C- 1 = -A, C^C- 1 = i$, (2.14) 

with a certain constant matrix C. This leads to a Z 4 automorphism in the present 
Hitchin system. This Z4 automorphism is not inherent in general Hitchin systems 
but is peculiar to the one describing the classical strings in AdS^. The system also 
exhibits Z 2 automorphism that corresponds to the reality of the string solutions. 



As argued above, n-cusp solutions are characterized by the condition that a — Y 
for | z | — ¥ oo. This is equivalent to the statement that one can always diagonalize $ 
and A at large \z\ by a suitable gauge transformation into the following form 

/ P(z) 1/4 \ 

, (2.15) 



hT l <S> z h -> -= 
y/2 



iP(zY' A 



\ 



-P(z) 1 / 4 



iP(z) 1 / 4 J 



h-'A z h + h-'dh -► m < ° 3 ° 



2 V o a 3 r 



h^Agh + h^dh-^-' 



m I a 3 



z \ a 3 



(2.16) 



where m and m are constants. 



2.3. Auxiliary linear problem and small solutions 

The system under consideration is classically integrable. This means that the con- 
nections appeared in the linear equations can be promoted to a set of one-parameter 
family of flat connections. That is to say, one can write down an auxiliary linear 
problem 



V z q(z,z-() = 0, Vjjg(*,Z;C) = 0, 



with 



Vi = D z + C 1 ^ 



D- z + ($- z 



(2.17) 



(2.18) 



where the Hitchin equations (12.121) . (12.131) are obtained as the compatibility condition 

[V«,Vi-] = 0. (2.19) 

One can estimate the asymptotic form of the solutions q(z,z;() for large \z\. 
There are four independent solutions, whose asymptotic forms are respectively given 
by 



e\z m z m ev^ 



(C -1 tu+Cu)) 



e 2 z- m z^e-^~ lw - C€,) , 



e 3 z m z m e V2 



(C _1 ^+C">) 



e 4 z m z m e^ 



((^w-Cw) 



(2.20) 



where w = j z P{z) l ^dz ~ z n ^ A and e*j are constant vectors. One can determine 
which of the four solutions shows the fastest decay for \z\ — > oo in each sector 



TTA 7r(2&-3) 4 A 7r(2/c — 1) 4 

Wk '■ 1 — arg Q < arg z < 1 — arg Q. 

n n n n 



(2.21) 



We call such solution the small solution Sk(z, z; () in the sector W k . Small solutions 
form a set of redundant basis of the solutions to the equations (j2.17p . We normalize 
the solutions so that 

(sj, Sj+i, s j+2 , s j+3 ) = 1, (2.22) 

where (si,Sj,Sk,si) = det(siSjS k si). As the differential operators are regular every- 
where in \z\ < oo, the sectors Wj and Wj +n are identified. Hence 

Sj-L-fi OC Sj. {Z.ZoJ 

One can determine the proportionality coefficient case by case whether n = 2k — 1, 
n = Ak — 2 or n = 4k for k G Z >0 . In the case of n — 6, 

s j+6 = ti^Sj (2.24) 

where 

ix = ± e 2 ™( m +™\ (2.25) 

Note that /x is a pure phase factor for solutions in the usual (3, 1) signature. For 
later use, let us also introduce the notation 

fj, = e* 4, (2.26) 

with being a real parameter. 

2.4. Y-system 

Making use of integrability, one can compute conserved quantities without knowing 
the explicit form of solutions. Instead, the fundamental building blocks are the Stokes 
data (si,Sj,Sk,si)((). Let us first recall some important identities concerning them. 
In addition to the normalization condition (I2.22p . there hold the following identities 

{s k ,s k+1 ,Sj,s j+ i)(C) = (sk-i } s k} Sj-i } Sj){iC) } (2.27) 

(sj, s k , s fe+ i, s k+2 )(0 = (sj, Sj_i, Sj_ 2 , s k )(i(), (2.28) 

which follow from the Z 4 automorphism. 

In [5] Alday, Maldacena, Sever and Vieira formulated the T-system and the Y- 
system associated to the general n-cusp solutions in AdS§. There, some particular 
Stokes data are chosen as T-functions and subsequently Y-functions are defined as 

6 



ratios of the form Y = TT/TT. Similar formulation of T- and Y-system was ob- 
tained from the spectral theory of ordinary differential equations [20]. By restricting 
ourselves to the n = 6 case, which is the simplest nontrivial case in their formulation, 
general formulas partly get simplified due to the relation (I2.22p . Y-functions are then 
defined bja 

Y x {9) = -( S2 , S 3,s 5 ,s 6 >(e e ), (2-29) 

Y 2 {9) = (s u s 2 ,s 3 ,s 5 )(s 2 ,s 4 ,s 5 ,s 6 )(e e+m / /k ), (2.30) 

Y 3 (6) = Y 1 (9). (2.31) 

These Y-functions are not entirely independent, but satisfy functional relations. By 
using Hirota bilinear identities (or Pliicker relations) among determinants and (12.241) 
(12371) . (T22HD, one can show that 

Yi ( 9 + t) Yi ( e - t) = l + Y2{el (2 ' 32) 

Y 4 + t) Y2 ( 9 ~t) =( 1 + A*nW)(l + Ai" 1 l r 1 W). (2-33) 

Another important property of the Y-functions is the periodicity 

Ya ( e+ ir) =Ya{9)) (2,34) 

which also follows from fF2T24"]) . fl2T27D . fT2T28|) . They also satisfy 



Y a {9) = Y a (-9), (2.35) 

which follows from the reality of string solutions. 

In fact, the Y-system obtained here is identified with that [TT] of Zjv-symmetric 
integrable models [131415] with N = A. Y a correspond to the fundamental representa- 
tions labeled by a = 1, 2, 3 of the A 3 Lie algebra. The present Y-system corresponds 
to the model with a chemical potential fi turned on. Note that the periodicity (I2.34p 
holds in the presence of \x. 

The cross-ratios of gluon momenta are given by special values of the Y-functions, 

^ Yx (^a\ Ut ^ +Y2 ((it±iM), (2 . 36 ) 



1 Y-functions introduced here are identified with those in [S] as Y\(6) = fi 1 [Y 1 A 1 MSV (ie 6 ')] 1 , 
Y 2 (9) = [Y 2 A ^(ie 6 )}- 1 , Y 3 (e) = n[Y 3 A ^(ie e )}- 1 . 



for k = 1,2, 3, where 

2 2 2 2 2 2 

C/i = hb 3 = X -lff, U 2 = b-A = X -f°f, U 3 = hb 2 = X -?ff. (2.37) 

x 13 x 46 x 24 x 15 x 35 x 26 

From the Y-system relations (J2.32J) . (12.331) one can verify that b/s are not entirely 
independent but obey the constraint 

bib 2 b 3 = 6i + b 2 + b 3 + [i + fi- 1 . (2.38) 

2.5. TBA equations 

The functional relations described above constrain the form of Y-functions mostly, 
but not entirely. To fully determine the Y-functions, we need an additional informa- 
tion such as the asymptotic behavior and the singularity structure. The asymptotic 
behavior of the Y-functions can be evaluated by the WKB analysis [21] ■ It can be 
shown that the Y-functions defined by (I2.29p - (l2.3ip exhibit the following asymptotics 



log 11(0) ->• \Z\e ±{ - e - l *\ logY 2 (0) ^V2\Z\e 



±{6-iv) 



for Re0^±oo, <p -~ < Im0 < ip + j, (2.39) 

where Z is related to the moduli parameter U in (I2.10p as 

Z ee \Z\e* = Ul J\l - e)Ut = ^^-Ul (2.40) 

As for the singularity structure, one can take Y-functions regular everywhere except 
at Re0 = ±oo. This is consistent with the functional relations (12 .32 p . (I2.33p . One 
could also consider Y-functions with singularities at finite 0. These correspond to 
excited states. Since we are interested in minimal area surfaces which correspond to 
the ground state, we restrict ourselves to the Y-functions regular at finite 0. 

Taking these into account, one can write down a set of integral equations which 
fully determine the form of Y-functions. By introducing the following notations 

e(9)=\ogY 1 (9 + ti P ), e(0) = logY 2 (0 + i<p), (2.41) 

the integral equations can be written in the form of TBA equations 

e = 2|Z|cosh0 + /C 2 *log(l + e~ £ ) + K x * log(l + /ie" £ ) (l + /^"V 6 ), (2.42) 

e = 2y/2\Z\ cosh + 2/Ci * log(l + e~ e ) + /C 2 * log(l + /ie" e ) (l + /i~V £ ), (2.43) 



where 



K,i{9) 



1 
27r cosh 9 ' 



Jc 2 (e) 



v^coshfl 
ir cosh 29 ' 



(2.44) 



and / * g = J d6'f(6 — 9')g{ff). Note that these equations are valid in the strip 
region — 7r/4 < lm.6 < tt/4. From these equations it is clear that e(9),e(9) are even, 
real functions 



-9) = e(9), 



-9) = ~e(9), 



e(6)=e(6), e(9) = e(9). 

In terms of Y-functions, these properties are expressed as 

Yai-O + ttp) = Y a (6 + i V ) 



(2.45) 
(2.46) 



(2.47) 



and (ESS]). 

Historically, the above TBA equations were first constructed for the Z4-symmetric 
integrable model. Indeed, the integral kernels can be obtained as derivatives of 
logarithm of the S-matrix elements in the Z/v model [13] with N = 4. 

By using f!2.36p . the cross-ratios (I2.37P are expressed in terms of e(9), e{9). If one 
wants to keep e, e evaluated within the above strip region, the appropriate choice of 
formulas are 



b k+l = exp 



kiri 



— up 



Uk-i = 1 + exp 



f(2k-l)m 

e| tip 



(2.48) 



for (2k - l)7r/4 < (p < kir/2, and 



7 fc+ i = exp 



kni 



up 



Uk = 1 + exp 



_/(2fc + l)7ri 



up 



(2.49) 



for kix/2 < ip < (2k + 1)tt/4, where k = 1,2,3 mod 3. The other cross-ratios are 
obtained by solving the relations ( I2.37P and f)2.38p . 

Solving the TBA equations numerically, we can calculate the cross-ratios as func- 
tions of \Z\,ip and <fi. In Figured] we plot 1/Uf. for various at fixed ip = — 7r/48 as 
an example o 
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-7r/48 is merely a generic value and does not have particular meaning. 




(a) 





(b) (c) 

Figure 1: (a) The \Z\ -dependence of the function 1/U\ for different values of at 
V = -tt/48. (b) l/U 2 . (c) I/I/3. 
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2.6. Area and remainder function 

As explained in J3j|5] in detail, the area of general n-cusp solutions with n 7^ 4k (k e 
^>o) 



Z>n) can be written a? 3 ' 



A = A div + A B DS-like + ^periods + -4free + COnst. (2.50) 

A div and A B DS-nke are contributions from the region inside the AdS radial coordinate 
cut-off and can be evaluated for general n. A per iods is given by a period integral over 
the Riemann surface of the ^-coordinate defined by dw = P{z) l ^dz. Af ree is shown 
to be identified with (minus) the free energy of the TBA system. 

On the other hand, there is known an all-order ansatz for the MHV gluon scatter- 
ing amplitudes proposed by Bern, Dixon and Smirnov [8]. It was shown (assuming 
the dual conformal symmetry) that BDS ansatz is correct for n = 4,5, but deviates 
from the string theory computation for n > 6 [9J. The function which complements 
the ansatz to produce the full gluon scattering amplitude is called the remainder 
function R. In the strong coupling limit, scalar magnitude of gluon scattering am- 
plitudes is then expressed as 

- A = -A div - A BDS + R. (2.51) 

A d [ v is identical to that in the string theory computation, with appropriate identifi- 
cation of the cut-off parameters. Abbs can be computed from one-loop perturbation 
and is known for general n. 

Taken altogether, the remainder function in the strong coupling limit is expressed 

as 

R = ^BDS — ^BDS-likc — ^periods — ^4frcc + COnst. 

= Ri- A P e r iods - A ivee + const., (2.52) 

where 

R\ = A BDS — A B DS-like- (2.53) 

Below we drop the constant term in (12.521) . 



For n — 4fc (k e Z >0 ) one needs extra terms [SJ[7]. 
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In the case of n = 6, terms in (I2.52p are explicitly given bye 

1 3 



fc=i 



.4 



periods 



(2.54) 
(2.55) 



-F = A frec = y j ^(2|^|cosh^log(l + /ie- e(e) )(l + ; u- 1 e- e(e) ) 

+ 2y/2\Z\ coshfllog(l + e~ m ) J . (2.56) 

Figure |5] and [3] show numerical results of the free energy and the remainder function at 
ip = — 7r/48 respectively, where the value of (p is the same as in Figured) From Figure 
[TJand[3Jb), we can read off the value of the remainder functions as the functions of 
the cross-ratios U^. 





(a) (b) 

Figure 2: (a) 3D plots of free energy A iree as a function of \Z\ and (0 < \Z\ < 

5, < < 37r/2) at ip = — 7r/48. (b) |Z|-y4f ree graphs for various values of at 
(p = -vr/48. 



Let us comment on choice of variables. General n-point gluon scattering ampli- 
tudes have 3n — 15 real moduli degrees of freedom. Correspondingly in the present 
case, the remainder function R for six-point amplitudes is a function of three inde- 
pendent real parameters. One can express R, either as a function of the momentum 
cross-ratios Uk or as a function of the moduli parameters \Z\,(p,<f>. The former choice 
respects the point of view of the four- dimensional M = 4 gauge theory while the lat- 
ter fit well with the two-dimensional description. In the following sections we mainly 



4 In 3 i?i is expressed as !?i = ^ fe=1 (| log Uu + \ Li2(l — '«/,)) with «/, = U k . This ean be 
rewritten as in (12.541) for Uk ^ (— oo, 0) by using a dilogarithm identity. 
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0.95 





(a) (b) 

Figure 3: (a) 3D plots of numerical data of remainder function R as a function of 

\Z\ and (f) at tp = — tt/48. (b) \Z\-R graphs for various values of at <p = — 7r/48. 



adopt the latter picture and study functional properties of R as well as its con- 
stituents A iree and R±, in particular in the two extreme regions \Z\ -C 1 and \Z\ 3> 1. 
The former picture in terms of the cross-ratios is also discussed. 

3. Free energy around CFT point 



The TBA equations describing the six-point gluon scattering amplitudes (I2.42p - 
(12.431) result from minimizing the free energy of the two-dimensional Z 4 -symmetric 
integrable field theory [TB] , which is obtained by deforming the Z4 parafermion CFT 
[T7] by the first energy operator [T4"l [To] . The model contains three particles with 
mass m, ylm and m, respectively. The third one is the anti-particle of the first. The 
free energy of the model on a space of length L ^> 1 with temperature 1/R gives the 
ground state energy of the model on a space of length R [12] . The TBA equations 
( E3Z]) - (g33D tel1 us tnat tne scale R is related to \Z\ as mR = 2\Z\. In [3], the free 
energy and the regularized area in the CFT limit mR — > are analyzed. In this 
limit, the cross-ratios take the values on the locus U\ — Ui — U3. In this section, we 
discuss the corrections to the free energy around the CFT point. 

We begin by noting that the action of the integrable model takes the form 



S = Spp + A / d 2 xe(x). 



(3.1) 



Here, SW is the action of the Z4 parafermion CFT, and e is the canonically normal- 
ized first energy operator with conformal dimension D £ = D £ = 1/3. The coupling 
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constant A > is related to the mass m as 



3'- 8 



3 



(2n\f=[m^l[-)\ 7(g), (3-2) 

where 7(2;) = T(x)/T(l — x). 

The free energy around the CFT point is then given by the ground state energy 



E(R) 
H2: 



of the perturbed CFT on a cylinder of circumference R with small A in (13.11) 



F = RE(R)-R 2 B(X), (3.3) 

where E(R) is expanded by the connected CFT correlators as 

{-X) n (2<K\1(D e --l)n+2 
n=l 



n=l 

/n 
(v(oo) e(z nj z n )--- e(z u z±) 1/(0) \ TJiz^f^dzl 

\ I connected ■*■ ■*■ 

i=2 

In the above, Eq is the unperturbed ground state energy, V is the operator corre- 
sponding to the vacuum, and the correlators are evaluated on the complex plane 
after a conformal transformation from the cylinder. We have also set z\ — Z\ — 1. 
Substituting A in (13.21) . one finds that the ground state energy has an expansion in 
(mR) 5™. Due to the Z2 symmetry e — > —e, only the terms with even n remain in 
the expansion [18]. 

The second term in (13.31) subtracts the bulk contribution to E(R), so that the 
free energy per unit length vanishes at zero temperature R — > 00 as implied by (12.561) 
and the TBA equations. This term is evaluated as [TB] 

B(X) = -\m\ (3.5) 

Although the derivation in [18] is given for \l — 1, the asymptotics needed there may 
hold also for [i 7^ 1 to give the same result. We will confirm that (13. 5 p is in agreement 
with the numerical results for both [i = 1 and /i^ 1. 

Without the chemical potential, namely, for /1 = 1, the vacuum operator V is 
the identity, and thus the expansion is straightforward. To proceed in the case with 
general /i, we recall the connection between the Z4 parafermion model and the spin-| 
XXZ (XXZ1/2) model, the Hamiltonian of which is 

N 

Hxxz = Y, [ S ! S !+i + SPI+i + cos « S j S Ui] ■ ( 3 - 6 ) 

3=1 



We have rescaled the free energy as R 2 F/L — >• F in accord with Q2.56p . 
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In the continuum limit, the spin variables are represented by a free boson $(r, a) 

<9 CT $ + (-l) j CiCOs($/r), 
c 2 cos($/r) + (-iyc 3 ], (3.7) 

where Sf = Sj ± iS? , c a are constants, $ is the dual boson, and (r, a) are the two- 
dimensional coordinates, r is the compactification radius, i.e., $ ~ $ + 27rr, and is 



S] 


r\j 




2ttt 


S I 


r\J 


i-rl> 

e2 



related to the coupling constant by r = v/2(l — ^) [21] in the unit where the selfdual 
radius is r s< j = a/2. Since the Z4 parafermion theory is described by a free boson 
compactified on the S 1 /^ orbifold with rpp = a/2/3 [25], the continuum limit of 
the XXZx/2 model with apF = |vr gives the Z4 parafermions. 

In terms of the XXZx/ 2 model, the chemical potential /i = e % ^ in the TBA equa- 
tions is understood as the twist parameter of the boundary conditions [261428] . 

S N+1 = S x , S N+1 = e 3 J S 1 . (3.8) 

The bosonization formula (13 .7p means that these are equivalent to the winding con- 
dition on the dual boson [21E9], $(r, o + R) ~ $(r, a) — -^dj, where R is the length 
of the space. For the original boson $, this induces the shift of the momenta, 

S-4< n -SF> (neZ) ' (3 ' 9) 

One can check that the ground state energy is also changed to (see also [30J) 

n ( 2d) 2 \ , 

6it V 3n (n — a) J 

due to the momentum shift, and that the TBA result in [3] is indeed reproduced for 

a = a-p-p: 

E = Ej =- —h-JL. . (3.11) 

Iq=q pf 6R\ tc 2 J v ' 

The above discussion shows that the vacuum operator V for general \x is identified 
with the momentum shift operator, 



^ = e -*sf?« = e -Vjr#* > (3.12) 

when the parafermions are bosonized. The energy operator e is also given simply by 
the vertex operators of $, since it is in the untwisted sector of the S ,1 /Z 2 orbifold 
compactification. The conformal dimension then determines its bosonized form, 



fl+ e'V3* + fl_e- ! V^ (3.13) 
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where a± are some coefficients including cocycle factors. Their explicit forms can be 
found in [3T] but, here, we only note that the products a±a T are c-numbers. 

Now, we are ready to evaluate the expansion of the ground state energy E(R). 
The first non-trivial term in ( 13. 4 p comes from the two-point function. Since e is 
canonically normalized as (e(z)e(0) )\cj>=o — M~ 4//3 5 one has for general <\> 

(e{z 2 )e{z x )) = (e i v / II*(oo)e i v / f*( Z2 )e- i v / f*( 2l ) e - i v / i^(o)^) 

= \l-z 2 \-*\z 2 \& 4 ', (3.14) 

with z\ = 1. The integral of this term is evaluated in (13.4p by the formula, 

d 2 u \u\ 2a \l - u\ 2b = vr 7 (l + 0)7(1 + &h(-l - a - b). (3.15) 

Collecting the results so far, we find the free energy near the CFT point to be 

F = BE. + |Z|" - C, 7 (i + IXl - ±) |Z|I + 0(|Z|¥), (3.16) 

where 

i\ n 



3 2 



1 /^ 

n 



7(4)] '7(^)7(3) ~ 0.18461. (3.17) 



1 8 



At = % the coefficient of \Z\~3 diverges, which implies that the CFT perturbation 
breaks down there. The expansions to higher orders are straightforward. Figure 
H] and show Af ree = —F from the above CFT perturbation and the numerical 
computation. The results in ( I3.14p and (I3.16P may be continued to imaginary <fi. 
While real <fi is physical in the context of the XXZ spin-chain, so is imaginary 
in the context of thermodynamics. Note that imaginary <p corresponds to minimal 
surfaces in (2,2) signature |3J. 

4. Remainder function around CFT point 

The remainder function (12.521) essentially consists of two non-trivial functions, A{ ree 
and R\. We studied the structure of A{ iee for \Z\ -^ 1 in the last section. Here let us 
see how R\ behaves as a function of \Z\,(p,<f> in the region \Z\ ^C 1. 

We know the form of R\ as a function of Uk (I2.54p . Using (I2.36P , one can express 
it in terms of Y-functions as 



ft ' Elt f.Jfi^)). (4 „ 



fc=i 
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Figure 4: (a) Plots of (\Z\, Af ree ) for various values of <fi at ip = — 7r/48. (b) Plots of 
(|Z|,A free ) for < \Z\ < 0.2 at <j> = vr/4, y> = -vr/48. 




Figure 5: Plots of A{ ree for < < 7r at |Z| = 0.1 and ip = — 7r/48. The dashed line 
represents the graph calculated by the formula (I3.16p . 
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What we need is the information of how \Z\,(p,<f> enters in the function Y 2 {6). From 
(I2.41I) - (I2.43I) we see that 9 and tp appear in the Y-functions only through the combi- 
nation 9 + iip. The Y-functions are periodic in 9 (12.341) and also exhibit Z2 symme- 
tries (12.351) . (12.471) . Also they are regular everywhere except at Re 9 = ±00, as we 
mentioned in section 2.5. Combining all these properties together, we see that the 
Y-functions admit the following Laurent expansion [TJ] 



Y a (9) = Y^(\Z\, ( l ) ) + J2Yl n \\Zim n + t- n ), t = ei^\ (4.2) 

n=l 

where Y} n E R. By substituting this into the expression (14. ip . one obtains the 
following Fourier expansion for Ri 



R 1 = J2 Ri i a) (\ z \^) cosAm V- ( 4 - 3 ) 



t ?( m )/ 

m=0 

Note that in deriving this expression, there occur cancellations among powers of 
t n evaluated at 9 = (2k + l)iri/4, k = 1,2,3, whenever n is not a multiple of 3. 
By construction the coefficients R™ are real. Therefore R± has shown to be real- 
valued, even, periodic function in ip with periodicity it/2. Our numerical results 
indeed exhibit this periodicity. 

Let us now analyze the form of R\ for \Z\ <C 1. It is seen from ( I2.42p - (l2.43p that 
the Y-functions grow dramatically at large 9 but do not vary so much in the region 
— log(l/|Z|) <C 9 -C log(l/|Z|). In other words, Y-functions draw a plateau over the 
region, which is wide when \Z\ is small. In order for the Y-functions to be so, the 
coefficients Y" a in the expansion ( 14 .2p have to be sufficiently small and at most 



Y a (n) (|Z|,0)~|Z|t n + (higher order terms) . (4.4) 

The appearance of powers of |Z| 4 / 3 = \U\ is not totally unexpected, as we have 
already seen in the last section that the CFT result for A {ree is obtained in powers of 
\Z\ 3. This suggests us that the Y-functions, and thus their descendants Uk and Ri, 
may well be expanded in powers of \Z\ 4 ^ 3 = \U\. Below we assume that this is the 
case. Let us then express Y 2 as 

00 
Y 2 (9) = Y,Yt\v + MA)\Z\l n . (4.5) 

n=0 

The first coefficient is the value of Y 2 evaluated at Z = and is known [3] as 

Y 2 (0 W) = 1 + ^ 2 /3 + yU -2/3 = 1 + 2c os^. (4.6) 
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(14. 4 p tells us that only the terms with n = 0,1 in (14. 2 p contribute to the next 
coefficient Y 2 . Therefore Y" 2 , as a function of cp, is at most the sum of a term 
independent of tp and a term proportional to t + t _1 = 2cos[A((p + i6)/3]. As a 
matter of fact, the former term is not allowed by the Y-system relations f !2.32j> . 
(12.331) . Hence Y 2 has to be of the form 



-(I) 



Yf } (<p,$)=yW(<l>)co8 



4;{ip + i6) 



(4.7) 



This form is also confirmed by numerical computations. At present we do not know 
the analytic expression for the function yW(0), but it can be evaluated numerically 
as in Figure Similarly, Y" 2 depends on cp only through cos[4(<£> + i9)/3] and 
cos[8(y9 + i0)/3], but actually we do not need the precise form here. By using the 
above data and the Y-system relations (I2.32p - (l2.33p . one can determine the behavior 
of -Ri for |Z| <C 1. If we express R\ in the form 



Ri = Y,ft?\<PMZ\* n , 



(U 



n=0 



the first few coefficients are obtained as 

R? W) 

Rf\vA) 



\u 2 {i 



4/3 2 ), 



0, 

3(4/3 2 



I-rlog(4/3 2 ))^ /(1)(0)2) 



64/3 2 (4/3 2 - If 



-y 



(4.9) 
(4.10) 
(4.11) 



with /3 = cos(0/3). Note that these three coefficients are (^-independent. This is 
consistent with the argument below (14. 3 j) that Fourier modes cos(4n(^/3) with n not 
being a multiple of 3 cancel out in R 1 . The 92-dependence could start appearing at 
the order of \Z\ . The above result agrees with numerical computation with high ac- 
curacy. Figure [7] shows a comparison between the above perturbative approximation 
with numerical plots. 

Collecting the results so far, the remainder function is expanded around \Z\ = 
as 



R 



1 - 



20 2 



7T 



4 



Li 2 (l-4/3 2 



+ 

+o(\z\ 4 ) 



-C8 7 - + — 
3 ' \ 3 3tt 



7t 3"37 



3(4/3 2 - 1 + log(4/3 2 )) 
64/3 2 (4/3 2 - l) 2 



y (1) (<Pf 



\z\l 



(4.12) 
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Note that \Z\ 2 terms in A periods and in A iTee cancel each other. Figure M shows a 
comparison between the perturbative approximation with numerical plots. 
One can also express the cross-ratios £4 (k = 1, 2, 3) as 

U„ = if + yU(<f,) (coB 4 ^-< 2 3 <:+1 >' r ) |Z|! + O (|Z|I) . (4.13) 

These are inverted to express \Z\,(p,<f> as functions of £/&: 

/3 2 = cos 2 ^ = 1(E/! + U 2 + U 3 ), (4.14) 

4 y/3(t/ 2 - Us) ,.„, 

t ' m 3 (p= 2U 1 -U 2 -U 3 > (415) 

4 -2^ + gg + U 3 

\ Z \ a = q a w ,x 3—- ( 4 -!6) 

dyW(0) COS 3<y9 

With help of numerical fitting, yW is also evaluated, e.g., as (see Figure [6]) 

y (1) (0) « 5.47669 - O.4841710 2 + O.O1194710 2 

w 1.31367 + 2.61136 cos | + 1.55402 cos -^. (4.17) 

Substituting these into (I4.12p gives an analytic expansion of R in terms of Uk, which 
can be directly compared with weak coupling results. This expansion describes the 
behavior of R around the locus U\ = U 2 = U 3 . 

One can check that the Jacobian of the above change of variables is proportional 
to (y^) 2 \Z\3 sin(|0) and the transformation is one-to-one for \Z\ ^ 0, < cp < 37r/2 
and < < 37r/2. The range of ip comes from that of the phase of U as seen in f)2.40p , 
while the range of <j) corresponds to the region where minimal surfaces are in (1, 3) and 
usual (3, 1) signatures [3]. From (I4.13J) . one also finds that all (\Z\,ip,(p) near the CFT 
point correspond to real cross-ratios. In the above approximation, each of \Z\,tp, 
has the following geometrical meaning in the parameter space (Ui,U2,Us): (14.141) 
implies that constant spans a plane perpendicular to the locus U\ = U 2 = U3. 
specifies the distance between the plane and the origin. On this plane, ip parametrizes 
a circle around the center U\ = U 2 = U3 = A/3 2 with a radius of a/3/2 y^((f))\Z\^ 3 . 
Thus, in terms of the cross-ratios, the weak dependence of R on ip for small \Z\ 
observed above is translated into that on the rotation around the locus U\ = U 2 = U3. 

5. Remainder function in large mass region 

So far, we have studied the remainder function when the mass scale \Z\ is small. In 
this section, following [12] we consider the large mass region, i.e., \Z\ ^> 1, which 
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Figure 6: Plots of y^((/>) obtained by fitting the numerical data of Y 2 (9) in the 
region 0.01 < \Z\ < 0.2 at tp = — 7r/48. When we fit the numerical data of y^((j>) 
by the function a + b<p 2 + c0 4 , we find that the coefficients are given by a = 5.47669, 
b = -0.484171, c = 0.0119471. If we fit them by a + 6cos(0/3) + ccos(20/3), we 
obtain a = 1.31367, b = 2.61136, c = 1.55402. 
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Figure 8: Numerical plots of \Z\-R vs perturbative solution of R (a) for various values 
of (f) at ip = -7r/48 (b) in the region < \Z\ < 0.2 at 4> = and <p = -7r/48. 



corresponds to the kinematical region around collinear limits. In this region, we can 
solve the TBA equations approximately. Throughout this section, we consider the 
case that — 7r/4 < <p < 0. Other cases can be analyzed similarly by using f)2.48|) . 

Let us first consider the free energy. For large \Z\, the pseudo-energies, e, e, 
behave as e(9) fh 2\Z\ cosh^ and e(9) ~ 2\/2\Z\ cosh#. The convolution terms in the 
TBA equations are suppressed exponentially. Thus the free energy is evaluated by 



A 



free 



d9 

2^ 



(ft + fi- 1 ) 2\Z\cosh9 e- 2 W coshe + 2V2\Z\cosh9 e- 2 ^ zlcosht 



(5-1) 



One can easily find that the free energy is expressed in terms of the modified Bessel 
function of the second kind, 



.4 



2\Z\ 



free 



71 



(fi + fi- 1 ) #i(2|Z|) + V2K 1 (2V2\Z\) 



(5.2) 



From the asymptotics of the modified Bessel function we see that A iree decays expo- 
nentially as \Z\ goes to oo. 

Next, let us consider the large \Z\ behavior of R\. The leading correction of b\ is 
given by 



^ _ e e(-*v) ~ e 2|2|cos(<p-7r/4)/j 



Si), 



(5.3) 
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where ip = ip + 7r/4 and 

(5.4 
The kernels are defined in (12.441) . Similarly, 

U 2 -l = e" { -^ « e 2 ^ z l cos ^(l + <5 2 ) (5.5 

with 

/oo ^ 

d0 2/C 1 ^ + ^)e" 2v/2|z|coshe +( ; u + /i- 1 )/C 2 ^ + 0)e" 2|z|cosh9 . (5.6 
-00 

From these and f)2.37|) . fj2.38|> . we can compute 62 and 63, 

63 = ^ = e^^+^l + 5,), b 2 = bl + I ? + fi + fi ~\ (5.7 

61 6163 - 1 

where 

<5 3 = — -= • (5.8 

1 + 01 

Thus the other cross-ratios U\ and U3 are given by 

1 - Ui = 1 - 6 2 6 3 = - e - 2v/2 l z l sin ^(l + Ai), (5.9 

l-U 3 = l- bxh = - e ^\z\*™$(i + a 3 ), (5.10 

with 

(1 + 5 3 + /ie - 2 l Z l cos ^+^ 4 ))(l + 5 3 + u-l e -2|2|co8(# +7 r/4))) 

Ai = — — 1, (5.11 

1 + C>2 
__ (1 + <$! + yUe - 2 l Z l cos ^-^ 4 ))(l + ^ + 0-l e -2|2|cos(0-V4))) 

A 3 = — — 1. (5.12 

1 + o 2 

Note that all 8i and Aj decay exponentially as \Z\ goes to 00. 

Here let us give a comment on the large \Z\ behaviors of the cross-ratios. When ip 
is far from zero, the cross-ratios show the asymptotic behavior (u\, u 2 , U3) — >■ (1, 0, 
as I Z I — > 00 where «& = 1/C4 (A; = 1, 2, 3). However if (p approaches zero, (ui, u 2 , M3 
can reach an arbitrary point on the segment (1 — c, 0, c) (0 < c < 1/2). To see this 
we need to take the double scaling limit (p — > +0, \Z\ — > 00 with \Z\ simp = o(> 
held fixed. In this limit the cross-ratios go to the point (see Figure [9]) 

(«!, u 2 , u 3 ) -> (1/(1 + e" 2 ^), 0, 1/(1 + e 2 ^)). (5.13 
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Figure 9: We take a large value \Z\ = 9.0 and plot the value of itj.( = 1/C4) for 
various ip at fixed = 0. Dashed lines correspond to the formula (15. 13j) obtained by 
taking the collinear limit. 

In order to invert the relation between (\Z\,(p,<j)) and the cross-ratios for large \Z\, 
one has to evaluate the integrals in (15. 4p and (15. 6p in more detail. 

To analyze the large \Z\ behavior of Ri, the following asymptotic expansion is 
useful: 

Li 2 (-e x (l + 5)) = -—- — - x5 + e~ x + 0{5 2 , e~ x 5, e~ 2x ) (x > 1 and S < 1). 

2 6 



(5.14) 



Therefore the asymptotic behavior of R\ is given by 

r 2 



R^ 



TT 



z + — 

1 ' 12 

+j(2y/2\Z\ cos<£ 6 2 + 2^2\Z\ sin^ A 3 - e - 2 ^l z l«»* + e -2y/2\z\On$ Ai y 

(5.15) 

The first term is divergent in the large \Z\ limit, but this term is canceled by the 
second term in (j2.52[) . Combining all the above results, we finally arrive at the large 
mass behavior of the remainder function 

7T 2 2\Z\ 



R 



12 71 



(/i + /I" 1 ) Ki(2|Z|) + v / 2i^i(2v / 2|^|) 



(5.16) 



+ ~(2>/2|Z| cos<£ 5 2 + 2^2|Z| sin^ A 3 - e -2^|coe* + e -i^\z\^ ^y 
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Figure 10: (a) Large \Z\ behavior of the remainder function for various <fi at <p — 
57r/24. Dashed lines show the asymptotic expansion obtained by (15.21) and (j5.15p . 
Points show the numerical results, (b) The |Z|-dependence of the remainder function 
vs the small (dashed line) and the large (dotted line) \Z\ expansions at <p = and 
if = — ir/ 48. The two expansions cover the whole region in \Z\ except a small region 
around \Z\ = 1.0. 



where 8 2 , A x and A 3 are given by (15.61) . (15.111) and (15.121) respectively. Since the 
second and the third terms decay exponentially to zero in the large mass limit, the 
remainder function approaches the constant 7r 2 /12. 

Figure [TD] shows the behavior of the remainder function in the large \Z\ limit. 
For the large \Z\, the asymptotic expansion of the remainder function (15.161) is in 
good agreement with the numerical data. Figure [TT] shows the (^-dependence of the 
remainder function for the large \Z\. We observe that (^-dependence is weak. This 
can also be seen in the analytic form (I5.16P where the (,3-dependence is suppressed 
exponentially in \Z\. 

6. Conclusions 

In this paper we have studied the remainder function of 6-point gluon scattering 
amplitudes in M = 4 super Yang-Mills theory by analyzing the TBA equations 
perturbatively and numerically. We have examined perturbative solution near the 
UV and IR limits and found that results are consistent with the numerical results. 
The remainder function is made of the free energy of the Z 4 -symmetric integrable 
model and the difference between the BDS part and the BDS-like part. The free 
energy near the CFT point can be obtained by the correlation functions with the 
chemical potential background. It is an interesting problem to generalize this result 
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Figure 11: ^-dependence of the remainder function (a) at \Z\ = 2.0 and = and 
(b) at \Z\ = 5.0 and = 0. The dashed line and points correspond to the asymptotic 
expansion and the numerical results, respectively. We see that R does not depend 
very much on (p. (R varies only up to 0.6% at \Z\ = 2.0 and 0.05% at \Z\ = 5.0.) 



to the case of n-point amplitudes corresponding to the generalized parafermions. 
The other part of the remainder function Ri is written as the sum of the dilogarithm 
function including Y-functions as arguments. In the present work, we could not 
completely determine its analytic form near the CFT point, but the undermined 
function of the chemical potential has been evaluated by numerical fitting. It is also 
an interesting problem to determine the series expansions of the Y-functions in order 
to know analytical properties of the remainder functions. In the large mass limit, 
we have obtained the first order correction to the remainder function. It would be 
useful to apply nonlinear integral equation approach [32] to analyze the TBA system 
further. 

Our results provide an analytic form of the remainder function for the 6-point 
amplitude at strong coupling away from the UV(CFT) and the IR(collinear) point. 
Such an analytic form beyond numerical ones will be important for further studying 
the super Yang-Mills theory at strong coupling. In particular, together with the 
recent results of the analytic form at weak coupling [33|IM]. which is still under 
active investigation, our strong coupling result will give a clue to understand the 
scattering amplitude to all order in the context of the AdS/CFT correspondence. 
In this regard, comparisons with the perturbative (both analytic and numerical) 
computations [TUJI33H35] would be of interest. One can expect that the physical 
picture of the amplitude to all order will emerge through further investigations both 
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at weak and strong coupling. As shown in Figure [TDJ we also find that simple first 
order expansions away from the UV and IR points give a good approximation to 
the remainder function for all the scale \Z\ as long as the expansions are valid. 
In addition, our results demonstrate that the identification of the two-dimensional 
integrable models underlying the four-dimensional super Yang-Mills theory [31 E] 
is useful, as well as interesting, for actual computations of the amplitude at strong 
coupling. The discussion in this paper may also be generalized to other TBA systems 
with chemical potential. 

The free energy Af ree is independent of (p, the phase of Z. From perturbative and 
numerical analysis, we observe that the 92-dependence of the remainder function is 
also weak: it starts appearing possibly at the order of \Z\ A for small \Z\ and at the 
order of e~ c ' z ' with c being some positive constant for large \Z\. In terms of the cross- 
ratios Uki the weak dependence for small \Z\ is translated into that on the rotation 
around the locus U\ = U 2 = U 3 . It would be important to further investigate the 
meaning of the 99- dependence to explore quantum corrections to the present analysis 
as noted in [5]. 
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